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Abstract 

In the present work, the notion of Cubic Spline Super Fractal Interpolation Function 
(SFIF) is introduced to simulate an object that depicts one structure embedded into an- 
other and its approximation properties are investigated. It is shown that, for an equidis- 
tant partition points of [xo,xn], the interpolating Cubic Spline SFIF g a (x) = g^\x) 
and their derivatives g&\x) converge respectively to the data generating function 
y(x) = y(°\x) and its derivatives y^\x) at the rate of h 2 ~^ +e (0 < e < l),j = 0, 1,2, 
as the norm h of the partition of [xq,xn] approaches zero. The convergence results 
for Cubic Spline SFIF found here show that any desired accuracy can be achieved in 
the approximation of a regular data generating function and its derivatives by a Cubic 
Spline SFIF and its corresponding derivatives. 
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1 Introduction 



Barnsley [T] introduced Fractal Interpolation Function (FIF) using the theory of Iterated 
Function System (IFS). Later, Barnsley [21 [31 H] introduced the class of super fractal sets 
constructed by using multiple IFSs to simulate such objects. Massopust [5] constructed super 
fractal functions and V- variable fractal functions by joining pieces of fractal functions which 
are attractor of finite family of IFss. 

FIF, constructed as attractor of a single Iterated Function System (IFS) by virtue of self- 
similarity alone, is not rich enough to describe an object found in nature or output of a 
certain scientific experiment. The objects of nature generally reveal one or more structures 
embedded in to another. Similarly, the outcomes of several scientific experiments exhibit 
randomness and variation at various stages. Therefore, more than one IFSs are needed to 
model such objects. A solution of fractal interpolation problem based on several IFS to model 
such objects is introduced in [6] by introducing the notion of Super Fractal Interpolation 
Function (SFIF). The construction of SFIF use more than one IFS wherein, at each level 
of iteration, an IFS is chosen from a pool of several IFS. This approach ensured desired 
randomness and variability needed to facilitate better geometrical modeling of objects found 
in nature and results of certain scientific experiments. 

Spline functions, introduced by Schoenberg [7], find vast applications in areas like data 
fitting [8], computer aided geometric design [HI [10], numerical solutions of differential equa- 
tions [IT], etc. A piecewise polynomial function d on an interval [xq, x^}, which is composed 
of subintervals aJj], i = 1, 2, . . . , N, is called a Spline of order n if (i) i)(x) is a continu- 
ous polynomial of degree atmost n — 1 in each subintervals [xt_i, Xi], % — 1, 2, . . . , N, and (ii) 
the derivatives $( m \ < m < n — 2, are continuous on [xq, Xn]. A Cubic Spline is a Spline of 
degree 3. For a data set {re,,} of n + 1 points, a Cubic Spline is constructed with n piecewise 
cubic polynomials between the data points. If $ represents a Cubic Spline approximating 
the function y e C 4 [x ,xn], then d is twice continuously differentiable and = y(xi). 
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Navascues and Sebastian [12] considered a Cubic Spline FIF as a generalization of classical 
Spline and obtained estimates on error in approximation of the data generating function by 
a Cubic Spline FIF. However, their Cubic Spline FIF was constructed using a single IFS and 
so it is not equipped enough to simulate an object that depicts one structure embedded into 
another. To approximate an object by a spline-like FIF, the concept of Cubic Spline SFIF 
is introduced in the present work and its approximation properties are investigated. The 
convergence results for Cubic Spline SFIF found here show that any desired accuracy can 
be achieved in the approximation of a regular data generating function and its derivatives 
by a Cubic Spline SFIF and its corresponding derivatives. 

The organization of the present chapter is as follows: In Section [2], a brief review on the 
construction of Super Fractal Interpolation Function for a given finite set of data is given. 
The notion of Cubic Spline SFIF is introduced in Section [3j It is proved in this section 
that, for an equidistant partition points of [xo,xjv], the interpolating Cubic Spline SFIF 
g a (x) = gf\x) and their derivatives g&\x) converge respectively to the data generating 
function y{x) = y(°\x) and its derivatives y^\x) at the rate of /i 2 ~ 3 ' +e (0 < e < 1), j = 0, 1, 2, 
as the norm h of the partition of [£ ,a;yv] approaches zero. 

2 Construction of SFIF 

In this section, a brief introduction on the construction of Super Fractal Interpolation Func- 
tion (SFIF) is given. 

Let So = {(xi,yi) G 1R 2 : i = 0,1,... TV} be the set of given interpolation data. The 
contractive homeomorphisms L n : / — > I n for n = 1, . . . , TV, are defined by 

L n (x) = a n x + b n (2.1) 

where, a n = and b n = x " x ^j x ^ _ For fe = 1, 2, . . . , M, M > 1 and n = 1,2, ... ,7V, 
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let the functions G n>k :Jxl->M defined by 

G n ,k(x, y) = in,ky + q n ,k( x ) ( 2 - 2 ) 

satisfy the join-up conditions 

G n ,k( x o,yo) = Un-i and G ntk (x N , y N ) = y n . (2.3) 

Here, 7^ are free parameters chosen such that yy n ,k\ < 1 an d ln,k 7^ ln,i for k ^ I. 

The Super Iterated Function System (SIFS) that is needed to construct SFIF corresponding 
to the set of given interpolation data So is defined as the pool of IFS 

{{R 2 ;u n , k :n = l,2,...,N}, = 1, 2, . . . , itf} (2.4) 

where, the functions u n<k '■ I x K — > I x R are given by 

w„,fc(x,y) = (L n (x),G n>k {x,y)) for all (x,y) e M 2 (2.5) 

By (12.3jl . it is observed that w ni fc are continuous functions. 

To introduce a SFIF associated with SIFS (J23}, let {W fc : i/(M 2 ) F(M 2 ), A; = 1,2,..., M}, 
be a collection of continuous functions defined by 

N 

W k {G) = [J u n>k {G), where u n ^ k {G) = {u n ^ k {x,y) for all (x,y) G G}. (2.6) 

n=l 

Then, {"H(M 2 ); Wi, . . . , W M } is a hyperbolic IFS, since h(W k {A),W k {B)) < max j nk h(A, B), 

l<n<N ' 

where h is Hausdorff metric on "H(IR 2 ). Hence, by Banach fixed point theorem, there exists 
an attractor A E n(n(M. 2 )). 
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Let A be the code space on M natural numbers 1, 2, . . . , M . In the construction of SFIF, 
for a a = 0\02 ... G A, let the action of SIFS (12.4p at the iteration level j be defined by 
Sj = W (Tj (Sj-.x), where So is the set of given interpolation data. For a fixed a G A, define 

G a = lim W Ck o . . . o W n (S ) = lim S k . (2.7) 

fc— >oo k— >oo 

The following proposition is instrumental for precise definition of a SFIF: 

Proposition 2.1 Let G a be defined by ( )2~71) . Then, G a is the attractor of SIFS (EH) for 
a = o\02 ■ ■ ■ Cfc • • • G A and is graph of a continuous function g a : I — > M. such that g a (x n ) = 
y n for all n = 0, 1, . . . , N. 

Super Fractal Interpolation Function (SFIF) is defined using Proposition 12.11 as follows: 

Definition 2.1 The Super Fractal Interpolation Function (SFIF) for the given 
interpolation data {(xi, yi) : i — 0, 1, . . . , A^} is defined as the function g a whose graph G a is 
the attractor of SIFS flU} . 

3 Cubic Spline SFIF 

The convergence of Cubic Spline SFIF is investigated here using the conditions of differ- 
entiability found in [6]. Throughout in this section, the interpolation data {(x n ,y n ) : n = 
0,1,..., A^} is assumed to be such that x n — x n _\ — h, for n = 1,2, . . . , N and x = 0. 
Also, throughout in the sequel, the SIFS (12. 4 p is chosen such that 7„ t fc = jk, n = 1,2, . . . , N, 
for some jk, < 7^ < 1, k G {1, 2, ... , M}. 

Definition 3.1 A SFIF g a , associated with SIFS ( |2T4|) . is called Cubic Spline SFIF if 

q n ,k, given in (12.21) . are cubic polynomials for n = 1,2, . . . , N and k = 1,2,..., M. 
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It is observed that, if q n , k (x) = q n ,k,3X 3 + q n ,k,2X 2 + q n ,k,\X + q n ,k,o, the coefficients q n)kji , i = 
0, 1, 2, 3, depend upon j k due to ( \2.3L necessitating in the sequel, the use of notation q n (lk, x) 
in place of q ntk {x). Throughout in this section, it is assumed that, for some A > 0, the 
polynomials q n satisfy 

\q n (7k,x)-q n (lfi,x)l <Ao / 31 x 



iTfe - 111 

for n = 1, 2, . . . , N, k, I — 1, 2, . . . , M and x £ [xo, Xn]. 

Let 7fc n be such that KfcJ < 1 , B kn = max I71 — j kn \ < 1 and q = ... q ... £ A be such 

1<1<M 

that qj = ko for all j £ N. Consider the family of continuous functions 

Q = {/ : / — > R such that / is continuous, /(x ) = f/o an d f( x N) = Z/Ar} (3.2) 



with metric dg(f,g) = max — For Read-Bajraktarevic operator T : A x Q — > Q 

defined by 

T(a,g)(x) = ffin \^G ik ^ k (l£(x), G ifc _i,«r fc -i o L£(x), G ik _ 2 , ak _ 2 ( 



G^L' 1 o . . . o L7\ x ),g(L- 1 0...0 L7\x))) ...)))} 



(3.3) 



where, A is the code space on M natural numbers 1,2, ... ,M and £7 is given by (13. 2p . the 
following proposition gives a bound on ||T(cr, <7)(#) — T(<^, 5»)(x)|| for a, q £ A: 

Proposition 3.1 Lei g E Q and inequality (13. ip fre satisfied. Then, 

\\T{a,g)-T{q, g )\\ oo <l3 k0 {y| 00 + ^ + _^-l (3.4) 
where a, q £ A, max |? re (7fc , x)| = -B , 7*= max |7/| < 1 and (3 ko = max |7«-7fc | < 1. 

xe^o.ajjv] 1<1<M 1<1<M 

n=l,2,...,JV 
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Proof By the definition of T(a, g)(x) (c.f. Q), 



|T( ff ^)(x)-r( ? ,^)(x)| 



k 



k / k \ 

+ E II W |^(7. P ,^ 1 o...o^r>))-?^(7^^ 1 °...°^r>)) 
p=i \j=p+i / 
fe-i / A; \ "i 

+ E II ^-Tfcol I?*,, (7*0, ^ o.-.oL- 1 ^)) I . (3.5) 
p =i \j= P +i J ) 

Since q n (lk , x) are cubic polynomials defined on compact set [x , xjv], there exists a B > 

such that max \q n (lk , x)\ = B . Therefore, by (13. ip and (13.51) . it follows that 

xe[x ,x N ] 



\T(a,g)(x) - T(q,g)(x)\ < lim | |j] \ lrT . - Tfco 

k / k \ fc-1 / fe-1 \ 

+ s ( n w j ^ - 7fc i ^0 + j2[ n - 7^1 ) b 

p=i \i=p+i / p=i \i=p+i / 

Since /3t n = max I7/ — 7i- n | < 1 and 7* = max I7/I < 1, ( 13.41) follows from the above 

K 1<1<M 1 ' ' 01 '* 1<1<M ]HI K V 

inequality. ■ 
The following proposition gives a bound on \\g a — g^W, a, q G A: 
Proposition 3.2 Let g a ,g f € and ( 13. ip be satisfied. Then, 

|| fc _ ft | U <_^_( WU + _A. + _*_) (3.6) 

where 7*,/?fc awd are as m Proposition \3.1[ 
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Proof Since (13. ip is satisfied, g a = T(a, g a ) and g q = T(q, g q ) for Read-Bajraktarevic oper- 
ator T, defined by (13. 3p . 

\\9a ~ &||oo < 7* \\g<r - &||oo + Pk ( ||&||o 



1-7* 1 - (3 ko 
The inequality (13. 6 p now follows from the above inequality. 

Remark 3.1 For 7fc = 0, 7* = /3& . Thus, inequality (13. 6 p for 7fe = implies 



1 I, . 7* / 4) + B 



1 - 7* V 1-7* 

By Hall and Meyer's theorem J73] /. ||g f ||oo — Koh 4 + Jo- Consequently, Proposition \3.2l for 
7fc = 9^es, 

\\g<r - svlloo < 7-^— (V 4 + Jo + • (3-7) 

1 - 7* V 1 ~ 7* / 

Using inequality (13. 7p . the order of approximation of data generating function y(x) by SFIF 
g a is given by the following theorem: 

Theorem 3.1 Let y(x) e C 4 [xq,xn] be a data generating function and g„ G Q be a SFIF 
associated with SIFS (12.41) such that %(h) = max |7j| < , for some s, < s < 1, 
where h = \xi — z = 1,2, ... ,N and \I\ = \x^ — Xq\. Then, for < e < s, 

\\y-g*\U = o{h 2 + e ). (3.8) 

provided (13.11) holds. 
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Proof Since (13. 1 p holds, an application of inequality (I3.7p gives, 



lb - 5UI00 < \\y - g f \\ + ||& - ^lloo 

i-7*(/i)V i-7.(/i)y 

1 /V ^4 , , M T , 7*W(4> + go) , 

-i^)r o/i +7 - WJo+ i-7,w (3 - 9) 

Using |7*(/i)| < the inequality (13. 9p implies 

II II ^ 1 J 1 (2+S) JV,4 , f JohM \ , AAq + Bq)^) -, 



The order of approximation error given by ( 13. 8ft follows from the above inequality. ■ 

Remark 3.2 It follows from inequality (13. lOj) £/iai ; in fact, \\y — g a \\oo — 0{h 2+s ). 

Remark 3.3 If M = 1 m STFiS 1 (12.41) . £/ien g CT reduces to a FIF. The convergence result for 
a Cubic Spline FIF fT^ follows as a particular case of Theorem \3.1\ 

The order in approximation of derivatives of data generating function by corresponding 
derivatives of SFIF is now investigated. It is known [6] that g& (x) and g£\x) are SFIFs 
associated with SIFSs |{M 2 ; uj ijk ,j(x,y) = G itk j(x, y)) : i = 1, 2, . . . , TV}, k = 

1,2, . . . , Mj for j = 1,2 respectively. Here, the functions Gi tkt i(x,y) and G i; fe )2 (a;,y) are 
given by 

Gi, k ,i(x, y) = N^ k y + Nq^i'jk, x) and G i>k , 2 {x, y) = N 2 j k y + N 2 qf\^ k , x). 



Let, for some Aj > 0, the polynomials q n satisfy 



|7fc - li\ 

for all n = 1, 2, . . . , N, k, I = 1, 2, . . . , M and x G [xq, xn\. 



<A 3 , j = 0,1, 2, (3.1i; 
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For j = 1,2, define the Read-Bajraktarevic operator Tj : A x Q — > Q by 



*fe-2i°"fe-2 J I *> • • • 



2i(<r, g)(x) = hm jCi^^j U^V), G ik _ ltffk _ ltj [l^ o L ik 1 {x),G i 

G iuauj (L7i 0...0... LT k \x), giL- 1 0...0... L~»))))) } (3.12) 

where, A is the code space on M natural numbers 1,2, ... ,M and Q is given by ( 13. 2p . To find 
the order of approximation of derivatives of data generating function y(x) by corresponding 
derivatives of SFIF the bounds on \\Tj(a, g)(x) — Tj(q, g)(x)\\ similar to ( 13. 4 p and the 
bounds on \\g& — g^\\,j = 1,2, similar to ( 13. 6 p are needed. Such a bound on \\Tj(a, g)(x) — 
Tj(q, g)(x)\\ for a, q G A, is given by the following proposition: 

Proposition 3.3 Let g £ Q and inequality (13. lip be satisfied. Then, for j = 1,2, 

WTfr^-TfrMKNifc |lkl|oo+ 1 /^ + 1 5^ fco | ( 3 - 13 ) 

where, a, q G A, max \q% \j ko , x)\ < B j} 7* = max \^\ < 

x(z[xo,x?{] 1<1<M 
n=l,2,...,N 

and (3 ko = max \j t - 7 fco | < ^. 
Proof By (l3~T2p . 



^ as, {^ fc (n ^ - ^oij 1^(4- 1 o . . . o l-, 

+ £ iv j(fc - m+1) f II It. 



2?) 



m=l \n=m+l 

X \q£ (l^LT^ 0...0 LT k \x)) - q% ( 7fco , L£o...o L"») \ 



k—1 / k \ 

+ E^ (fc " m+1) II \^n-lK \) \qi2(% ,L^o...oL-\x))\\. (3.14) 

m=l \n=m+l / ^ 
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Since qn\^k Q -,x) are polynomials denned on compact set [io,%], there exists a Bj > such 

that max (lk ,%)\ < Bj. Therefore, by (13. lip and (j3.14p . it follows that 

xe[x ,x N ] 

n=l,2,...,7V 



iTj^g^-T^g^x)] < lim lN jk f J] 1 7^-7*0 1 J IMU 

m=l \n=m+l / 

m=l \n=m+l / ^ 



Since /3fc„ = max (hi — 6\) < kfs and 7* = max I7/I < -rW, 03.131) follows from the above 

u l</<Af JV 1<1<M JV 

inequality. ■ 

,0) Jj) 1 



Similar to Proposition ^. 2\ the following proposition gives a bound on || go- — g\ || for cr, q G A: 
Proposition 3.4 Lei Ocr,^ G and ( 13. lift be satisfied. Then, for j = 1,2, 

|| 9 0) _ „ )| U < jgflj. J,, M | U + A, NIB, 1 

y? 1100 - l-i\^7* 1100 1-A^ 7 * l-Wp^ J v ; 

where, 7*,/3fc and £?,,• are as in Proposition ^. 

Proof Since (13 . 1 1 [) is satisfied, = Tj(a,ga) and g<? = Tj(q,gq) for Read-Bajraktarevic 
operator Tj, defined by ( 13.121) . 

W - 9®\U < Nj l*W$ ] ~ 9?\U + N j p k0 (ll^lU + + ; N3 £; R 

[ 1 - iv J 7* 1 - ^ J Pifc 

The inequality (13. 15ft now follows from the above inequality. 
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Remark 3.4 For 7fc = 0, 7* = /3& . Thus, inequality ( 13 . 1 5 j) implies 



^ 1 1 00 — 1 Ar^, 1 life H°° ^ 



5y i/a// and Meyer's theorem pSy . \\g¥ \\oo < Kjh + Jj. Consequently, Proposition \3.4\ for 

lk = 0, gives, 

ll^-^ll < ^ fl^ 4 I J I Vb^ l (3 16) 

iif CT <7 S iioc< ^^1^ + •>,+ j- ^- ib ; 

Using inequality ( I3.16p . the orders of approximation of derivatives of data generating function 
y(x) by corresponding derivatives of SFIF g a are given by the following theorem: 

Theorem 3.2 Let y(x) G C 4 [xq,xn] be a data generating function and g a G Q be a SFIF 
associated with SIFS (12.41) such that 7*(/i) = max |7j| < j^fj, for some s, < s < 1, 
where h = Xi — i — 1, 2, . . . , N and \I\ = x^ — xq. Then, for j = 1, 2 and < e < s, 

Wv^-gP Woo = o(h 2 -i + <) (3.17) 

provided ( 13. lip holds. 

Proof Since (13. lip holds, an application of inequality ( I3.16P gives, 

||„(i) _ „0')|| 
1 1 y Her 1 1 00 

<\\y i ' , -'A"\\ + U"-g { J , \\^ 

< -T^WW){ Klh " , + Nl '" {h) J ' + T^k to^'BA \ (3.18) 
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Using |7*(/i)| < | jja+g , the inequality (13. 18[) implies 



\\y {3) - g^Woc 

- |J|(2+8-i) _ " + |J|(2+»-i) + |J|(2+«-i)_/ l (2+«-J) L^ + iV 

(3.19) 

The order of approximation error given by (13.171) follows from the above inequality. ■ 

Remark 3.5 It follows from inequality (13.191) that, in fact, \\y^ — gz ■ ||oo — 0(h 2+s ~^), 
J = 1,2. 

Remark 3.6 J/M = 1 in S/FS |{R 2 ; w^ac, y) = (£;0), G^O, y)) : i = 1,2, ... ,N},k = 
1,2, ... , M j, j = 1,2, then g&\ j = 1,2, are FIFs. The convergence results for derivatives 
of a Cubic Spline FIFs IW^ follow as a particular case of Theorem \3.S[ 



4 Conclusions 

In the present work, the notion of Cubic Spline SFIF is introduced for an efficient approxima- 
tion of the data generating function. The approximation properties of a Cubic Spline SFIF 
are investigated and it is proved that the order of approximation of the data generating func- 
tion y(x) = y(°\x) and its derivatives y^\x) by an interpolating Cubic Spline SFIF g a = g^ 1 
and its corresponding derivatives gi respectively, is o{h 2 ~i +e ), < e < 1, j = 0, 1, 2. These 
convergence results show that it is possible to approximate any regular data generating func- 
tion by Cubic Spline SFIF with arbitrary accuracy. Our study of Cubic Spline SFIF is likely 
to have wide applications like pattern-forming alloy solidification in chemistry, blood ves- 
sel patterns in biology, signal processing, fragmentation of thin plates in engineering, stock 
markets in finance, wherein significant randomness and variability is observed in simulation 
of various processes. 
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